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We derive the pre-shock density and scale length along the line of sight for 
the collisionless shock from a deep HST image that resolves the Ha filament in 
SN1006 and updated model calculations. The very deep ACS high-resolution 
image of the Balmer line filament in the northwest (NW) quadrant shows that 
0.25 < n < 0.4 cm -3 and that the scale along the line of sight is about 
2 x 10 18 cm, while bright features within the filament correspond to ripples with 
radii of curvature less than 1/10 that size. The derived densities are within the 
broad range of earlier density estimates, and they agree well with the ionization 
time scale derived from the Chandra X-ray spectrum of a region just behind the 
optical filament. This provides a test for widely used models of the X-ray emission 
from SNR shocks. The scale and amplitude of the ripples are consistent with 
expectations for a shock propagating though interstellar gas with ~ 20% density 
fluctuations on parsec scales as expected from studies of interstellar turbulence. 
One bulge in the filament corresponds to a knot of ejecta overtaking the blast 
wave, however. The interaction results from the rapid deceleration of the blast 
wave as it encounters an interstellar cloud. 
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SN1006 (G327.6+14.6) is one of the best 
SNRs for studying the physics of colli- 
sionless astrophysical shocks, in particu- 
lar the acceleration of non-thermal parti- 
cles. It is a nearby Type la supernova 
remnant at a distance of 2.1 kpc (the dis- 
tance we assume throughout) with a diam- 
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eter of ~ 18 pc (IWinkler. Gupta k Long 



20031 ) and a shoc k speed in the 2500 



2900 km s 1 range ( 


Ghavamian et al. 


2002; 


Heng k McCrav 


2006). The remnant has 



a high Galactic latitude and modest fore- 
ground reddening, E(B-V)=0 11 ± 0.02 
(ISchweizer k Middleditchlll980h . 



nent profile (j Chevalier k Raymond! 1 19781 ). 
The broad component is due to charge 
exchange between neutrals and protons, 
which produces a population of neutrals 
at nearly the post-shock proton tempera- 
ture. The narrow component is produced 
when cold ambient neutrals pass through 
the shock and emit line radiation before 



SN1006 has been observed at radio 
(iRevnolds k Gilmore!ll993l:lMoffett. Goss k Revrijtlj ge transfer or ionization occurs. The 



1993 1), optical (Ghavamian et al 



2002 



ratio of the broad to narrow flux is sensi- 



Kirshner. Winkler fc Chevalierlll987l:lsmith et alltive to the electron and ion temperatures. 



1991 ). ultraviolet (IRavmond ] Blair k Lon 



199 51: IKorreck et al. I I2004T) and X-rav 



([Koyama et al!l995l:IWinkler. Gupta k Lone! 



200 31: ILong et a 



2003: Bamba et al. 



2003 



Dyer et al.ll2004j ) wavelengths. Pure Balmer 



line filaments w e re fou nd in the optical by 



van den Bergh I (119761 ). In the radio and 



X-ray, the remnant has a limb-brightened 
shell structure with cylindrical symme- 
try around a SE to NW axis probably 
aligned wit h the ambient g alactic mag- 
netic field ( iReynoldsl Il996l ; IJones k Pye 
1989h . The NE shock front of SN1006 



shows strong n o n-the r mal X-ray emission 
(kovama et ail 1 19951 : iDver et all l2004f ) 
while the NW shock shows very little 
non-thermal emission at radio or X-ray 
wavelengths. On the other hand, the pre- 
shock density is sev eral times higher in th e 
NW than the NE (IKorreck et al. 1 120041 ) . 
Knots of X-ray emission from shocked SN 
ejecta are sca ttered through the interior of 



the r emnant (ILong et al.l 120031 : IVink et al. 



2003|). 



A Balmer line filament defines the blast 
wave in the northwest quadrant of SN1006, 
and the Ha profile provides diagnostics for 
the shock speed and ion-electron thermal 
equilibration. The Ha emission from a 
Balmer-dominated shock has a two compo- 



The FWHM of Ha broad component is 
2290 ± 80 km s _1 , and models imply a 
shock speed of v shock = 2890 ± 100 km s _1 
for a shoc k with little electron-i on equi- 
libration (jGhavamian et al.l 120021 ) . How- 
ever, new models incorporating some addi- 
tional physics obta in a lower shock speed o f 
2509±111 km s" 1 (iHeng k McCravll2006h . 

The pre-shock density, n , is an im- 
portant parameter for understanding the 
evolution of SN1006, as well as for in- 
terpreting the X-ray spectra in terms 
of ionization time scale and determin- 
ing the relative contributions of shocked 
ISM and SN ejecta to the X-ray emis- 
sion. The density could also be im- 
portant for attempts to understand the 
high ratio of non-thermal to thermal X- 
ray emission in this SNR. Estimates of 
n cover a wide range, from 0.05-0.1 
cm" 3 based on the global X-ray emission 
( jHamilton. Sarazin k Szymkowiaki 119861 ) 
to 1 cm" 3 based on interpreting the scale 
over which the X-rays brighten as the 
length s cale for ionization of t he shocked 
plasma (IWinkler k Long||l997l ). A related 
parameter is the length scale along the 
line of sight. The ripples in the SNR 
blast wave could in principle result from 
density inhomogeneities in the ambient 
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medium or from knots of SN ejecta over- 
taking the blast wave. ISM density fluc- 
tuations have been inferred for a section 
of the non-radiative sho ck in the Cygnus 



Loop ([Raymond 1 120031 ) . while an ejecta 



knot is clearly the cause of o ne bulge in 



the Ha filament of S N1006 flLong et al. 



20031 : IVink et al.ll2003h . The issue is im- 



portant for estimating the amplitude of in- 
terstellar turbulence on sub-parsec scales 
and for the interpretation of the distance 
between the blast wave and the reverse 
shock in terms of particle acceleration 
flWarren et al.l l2005h . 

In this paper we use an Ha image ob- 
tained with the ACS imager on the Hub- 
ble Space Telescope to determine the pre- 
shock density and the length scale along 
the line of sight. This is possible because 
the ACS images resolve the thickness of 
the narrow zone behind the shock where 
hydrogen atoms are excited and ionized, 
and that thickness scales inversely as the 
pre-shock density. We have computed new 
models of the Ha emissivity as a function 
of distance behind the s hock taking into 
accou nt recent results by iHeng &: McCray 
(120061 ). We discuss the observations in 
the next section, then compare with mod- 
els to derive the shock parameters. Sec- 
tion 4 provides a limit on the brightness of 
any shock precursor, compares the densi- 
ties with other estimates and discusses the 
implications for other analyses of SN1006 
observations. Section 5 summarizes the 
conclusions. 

2. Observations 



The Hubble Space Telescope's ACS 
Wide Field Camera (WFC) imaged a full 
field of 202 x 202 arcsec 2 at coordinates 



a 200 o =15 h 2 m 19.02 s , 5 2000 =-41° 44' 48.4" 
on 9 orbits from 15-17 February 2006. Ex- 
posures with durations of 2,746.0 s, 2,848.0 
s, and 2,828.0 s were obtained for each sub- 
set of 3 orbits, for a combined 25,266 sec- 
ond exposure. The exposures were taken 
with the F658N H-alpha filter, which has a 
flat response at wavelengths above about 
6558A, drops to half the peak transmis- 
sion at about 6548 A and 20% of the peak 
transmission at 6540 A. This means that it 
passes the narrow component, all the red 
wing of the broad component, and about 
half the blue wing of the broad component, 
or about 90% of the Ha emission. 

The image was centered on the position 
where Ghavamian et al. (2002) obtained 
a low dispersion spectrum and Sollerman 
et al. (2003) obtained a high dispersion 
spectrum, so we are able to use the shock 
speed, neutral fraction and electron-ion 
equilibration derived from those observa- 
tions. The position was also observed in 
the ultraviolet with the Hopkins Ultra- 
violet Telescope ([Raymond. Blair fc Long 
19951: lLaming et al.lll996 h and with FUSE 



(IKorreck et al. 1120041 ). so that we can use 



the information derived in those papers 
about ion-ion thermal equilibration to con- 
strain model parameters. The Ha filament 
also defines the outer boundary of posi- 
tion NW-1 in the X-ray spectral analysis 
of Long et al. (2003), who found the spec- 
trum to be consistent with thermal emis- 
sion from shocked interstellar gas. 

The raw data were combined and re- 
duced with the standard ACS calibration 
pipeline (CALACS), involving bias/dark- 
current subtraction, flat-fielding, image 
combination, and cosmic-ray rejection 
(Sirianni et. al. 2005), then 'drizzled' 
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(Fruchter and Hook 2002) onto a 0.03 arc- 
sec pixel scale with the task 'multidrizzle' 
to correct for geometric distortion and im- 
prove th e sampling of the poin t spread 



function (IKoekemoer et al. II2002I ). 



Figure [T] shows the HST image overlaid 
on the Chandra X-ray image. The mor- 
phology of the filament is clearly that of a 
rippled sheet seen edge-on, with the bright 
rims corresp onding to t angencies to the 
line of sight (lHesterlll9871 ). For the present 
purposes we are interested in the simplest 
tangencies, since those are amenable to 
modeling. The bulge near the SW cor- 
ner of the image is morphologically sim- 
ilar to a larger, brighter region farther to 
the SW where a clump of ejecta is overtak 



ing th e shock ([Long et al.ll2003l ; IVink et al. 
20031 ). but it shows only a slight X-ray en- 



hancement. The shock morphology in the 
bulge is more complex, so this paper con- 
centrates on the smoother regions of the 
filament. 

IDL was used to extract and plot the 
curve of the shock front spanning the driz- 
zled image and to find the approximate 
direction perpendicular to the shock. At 
each of 40 positions we extracted the spa- 
tial profiles across the shock for a range 
of angles near the initial estimate and se- 
lected the the profile showing the narrow- 
est Ha peak as the one closest to the shock 
normal direction. Many of the profiles suf- 
fer from low signal to noise or from com- 
plexity due to several tangencies to the 
line of sight. We have selected 8 profiles 
with bright, simple Ha peaks for further 
analysis. In particular, we model sections 
of the trailing edge of the filament in re- 
gions corr esponding to sections F. G and 
H of the IWinkler. Gupta k Lond fl2003h 



proper motion analysis. Figure [2] shows 
the boxes used to extract spatial profiles, 
starting with profile 8 in the upper left 
and ending with profile 29 in the lower 
right. The profiles were extracted using 
boxes 4.5"wide and 9"long. Their positions 
and position angles are shown in Table 1. 
The pro files extracted were similar to th ose 
used by ( IWinkler. Gupta fc Longil2003l ) for 
measuring the proper motion of the fila- 
ment; here, however, the goal was to de- 
termine the width, not the position, of the 
filament. 

3. Analysis and Results 

To interpret the images we compute 
models of the Ha brightness behind a 
curved shock, convolve the model intensity 
distribution with the ACS point spread 
function, and compare the models to the 
observations. In general, the pre-shock 
density controls the brightness drop off be- 
hind the peak, because the thickness of the 
emission region is inversely proportional to 
the density. The radius of curvature of the 
shock determines the fall off ahead of the 
peak for the concave outward filaments 
that we model here. The absolute inten- 
sity scales approximately as nlfneutR 1 ^ 2 , 
where n is the pre-shock density, f neut is 
the pre-shock neutral fraction, and R is 
the radius of curvature of the shock. 

3.1. Model Calculations 

We start with model calculations simi- 



lar to those of lLaming et al.l (119961 ). Fig- 
ure [3] shows the Ha emissivity SIS db func- 
tion of distance behind the shock for a 
2900 km s 
density n 

tion of 0.1, and a ratio T e /T p 



1 shock model with pre-shock 
= 0.25 cm -3 , a neutral frac- 

= 0.05 
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at the shock front. The model follows 
neutral hydrogen as it passes through the 
shock and undergoes collisional excitation 
and ionization by protons and electrons as 
well as charge transfer with post-shock pro- 
tons. The at o mic r ates are described by 
Laming et al.l (119961 ) . and Coulomb colli- 
sions slowly transfer energy from the ions 
to the electrons. The radiative transfer in- 
volved in the conversion of Ly/3 photons 
to Ha photons, important f or the narrow 



compo nent, is described in lLaming et al. 
(119961 ) . A fraction of the Ha arises from 



converted Ly/3 photons, and those Ha pho- 
tons are produced over a scale of about 1 
Ly/3 mean free path. That scale is about 
0.1", which is small enough compared to 
the observed filament widths that we ig- 
nore it. 

Evidence for a shock precursor has been 
reported for a number of non-radiative 
shocks from low ionization emission lines 



( He ster. Raymond fc Blairlll994j ; ISollerman et al. 
20031 ). from the veloc ity widths o f their 
narrow components (Smith et al.l Il994l ; 
Hester. Raymond fc Blairlll994 ). and from 
the spatial distrib ution of the narro w com- 
ponent emission (ILee et al. 1 120061 ). The 



narrow component in SN1006 shows no 
broadening be yond that expected for the 
ambient ISM (ISollerman et al. 1 120031 ) . so 
we do not include any precursor emission 
in the models. 



Recently, iHeng fe McCravl (120061 ) inves- 
tigated the consequences of the sharp de- 
cline of the the charge transfer cross sec- 
tion at speeds above about 2000 km s _1 
for the Balmer line profiles. They esti- 
mated a shock speed of 2509±111 km s _1 
fro m the Balmer li n e pro file presented 
by iGhavamian et al.l (120021 ). as opposed 



to the value 2890±100 km s" 1 found by 
Ghavamian et al. Both values apply to 
the case of little ion-electron thermal equi- 
libration in the shock. Strong equilibra- 
tio n is ruled o ut by the X-ray spectrum 
of lLong et al.l (120031 ) and by the combi- 
nation of b road to narrow intensity ra- 
tio given by IGhavamian et al.l (120021 ) and 
the Ly/3 radia t ive tr ansfer calculations of 



Laming et al.l (119961). The smalle r shock 



speed of IHeng fc McCravl (120061 ) would 
decrease the distance scale of the Ha 
emission behind the shock. However, it 
would also imply a smaller distance for 
SN1006 based the combination of the 
proper motion and the and s hock speed 
(jWinkler. Gupta fc Longll2003l ). so the an- 
gular width of the emission region would 
be unchanged to first order. 

There is anot her important implica- 
tion of the work of IHeng fc McCravl (120061 ) 
for the present study. The rapid drop in 
charge transfer cross section with increas- 
ing velocity means that neutrals are more 
likely to undergo charge transfer with pro- 
tons moving away from the shock than 
with protons moving toward the shock. 
While the velocity dependence of the cross 
section was included in earlier models, the 
anisotropy of the resulting H I velocity dis- 
tribution was not. Thus the earlier models 
implicitly assumed that the broad compo- 
nent neutrals move away from the shock 
at Vs/4, the same speed as the post-shock 
protons. By integrating the product of ve- 
locity tmiej_charge transfer cross section, 
a cx v ( Barnettl 1990 ). over Maxwellian dis- 
tributions at the post shock temperatures, 
we find that after 1 charge transfer event 
the average neutral is moving away from 
the shock at 1500 km s _1 behind a 3000 
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km s -1 shock, rather than the 750 km s _1 
of the ionized gas. For a 2500 km s _1 
shock, the neutrals move at 1100 km s^ 1 
instead of 625 km s^ 1 . Thus for the inter- 
esting range of shock speeds, the neutrals 
move away from the shock twice as fast as 
is assumed in the model shown in Figure 
[3j The relative velocity is much smaller 
for the second charge transfer event, so 
after two charge transfers the neutrals 
have a speed closer to that of the down- 
stream plasma. Kevin Heng (2007, private 
communication) has provided the average 
downstream speed of the broad component 
neutra ls computed by the lHeng fc McCray 
fl2006h model code. For the 2500 to 3000 
km s _1 velocity range it is 1.32 times the 
plasma speed. Therefore, for comparison 
to observations we stretch the spatial scale 
of the broad component emission shown 
in Figure 3 by a factor of 1.32. Particle 
conservation implies that the neutral den- 
sity is decreased by the same factor, so the 
broad component emissivity at each point 
is reduced by a factor of 1.32. Finally, 
for comparison with the observations we 
multiply the broad component emission 
by 0.75 to account for the drop of sensi- 
tivity in the blue wing resulting from the 
transmission of the F658N filter. 

To model the geometry of the filament, 
we assume a shock surface that is concave 
outward with radius of curvature R. The 
Ha emissivity from the planar model ex- 
tends behind the shock at each point. Nu- 
merical integration then gives the bright- 
ness as a function of position relative to 
the tangent point of the shock. For com- 
parison with the observations, we convolve 
the model emissivity with the ACS point 
spread function. In order to properly ac- 



count for the pixel sampling in cuts made 
at about 45° to the rows and columns 
of the detector, we measured the Gaus- 
sian widths of stars near the filament (2.31 
pixel FWHM), then placed a dense series 
of circular Gaussians along a 45° line, ex- 
tracted the profile perpendicular to that 
line, and measured a slightly broader 2.43 
pixel width. This Gaussian is convolved 
with the models for comparison with the 
observed spatial profiles. 

Figure H] is a schematic diagram of the 
geometry we imagine in a plane that in- 
cludes the line of sight (dashed line) and a 
radial vector from the center of the SNR. 
The light line shows a large scale ripple in 
the shock front, and we have superposed a 
smaller scale ripple to obtain the shape of 
the heavier curve that is shaded toward the 
inside of SN1006. The shading indicates 
the Ha emissivity, which fades gradually 
from the shock front towards the inside of 
the SNR. The trace in the lower right indi- 
cates the Ha brightness obtained for a cut 
across the tangent point of the shock by in- 
tegrating the Ha emissivity along different 
lines of sight. 

3.2. Comparison of Models and Ob- 
servations 

Figure [5] compares a grid of models 
with different n and R to the spatial 
profile at position 10. The sharp spikes 
ahead of the Ha peak are faint stars in 
the extraction region. The models have 
been scaled to match the peak Ha bright- 
ness of the filament by simply adjusting 
the neutral fraction. Neutral fractions 
above 1 are obviously unphysical, so mod- 
els with low n n and sma ll R are ruled out. 
Ghavamian et al.l (120021 ) estimated a pre- 
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shock neutral fraction of 0.1 from the ratios 
of He I and He II lines to Ha, so we take 
the permitted range of neutral fractions to 
be 0.05 < f neu t < 0.2. Figure [5] shows 
that models with uq below about 0.30 cm -3 
fall off too slowly behind the shock, while 
the models with n above about 0.35 cm -3 
give too sharp a peak. The R=5 x 10 16 
cm model comes closest to matching the 
observed profile. Models with smaller R 
predict Ha emission fainter than observed. 
We conclude that 0.3 < n < 0.35 and 



10 



16.5 



< R < 10 



16.8 



cm. 



Figure [6] shows the analogous plots for 
position 28. The peaks in this section of 
the filament are both brighter and broader, 
with a fairly constant Ha intensity ahead 
of the brightness peak. None of the mod- 
els match exactly, probably because the 
rippled sheet does not follow the assumed 
shape of an arc of a circle. From Figure M 
we conclude that no must be greater than 
0.25 cm~ 3 to avoid a long tail toward the 
inside of the remnant, but that densities 
above 0.4 produce too sharp a peak. R 
must be less than 2 x 10 17 cm, because 
larger radii of curvature predict a shoul- 
der on the outer side of the spatial pro- 
file that exceeds the observations, while 
R less than 0.7 x 10 17 cm requires unac- 
ceptably high values of f neut to match the 
brightness. We conclude that the accept- 
able ranges are 0.25 < uq < 0.35 and 
0.7 xlO 17 < R < 1.5xl0 17 cm. We note, 
however, that the agreement between the 
model spatial profiles and the observations 
is not as good as at position 10. Based on 
Figure [2] , it seems possible that the bright 
filament at position 28 contains more than 
one tangency to the line of sight, broad- 
ening the spatial profile and increasing the 



total brightness. Thus the upper limit on 
n and the lower limit on R are less secure 
than in the case of position 10. 

Table 1 presents the ranges of n and R 
derived from other positions along the fila- 
ment. We did not attempt to model other 
positions where the emission peak is faint 
or where complex morphology indicates a 
more complex structure than can be ap- 
proximated by a simple curved sheet. The 
profiles of positions 8 through 11 are qual- 
itatively similar to that of position 10, and 
those of positions 26 through 29 are like 
that of position 28. 

4. Discussion 

4.1. Pre-shock density 

The densities near 0.3 cm" 3 that we de- 
rive fall in the middle of the range of previ- 
ous va lues. lHamilton. Sarazin fc Szymkowiak 
( 119861 ) estimated 0.05-0.1 cm -3 from the 
global X-ray spectrum. Since the pre- 
shock density in the NE, and quite likely 
the rest of the remnant, is about 2. 5 
(lLong et al.1 l2003h to 4 flKorreck et al 



2004 ) times smaller than in the region we 
observ ed, the lHamilton. Sarazin fc Szymkowiak 
( 119861 ) estimate is consistent with our re- 
sults, even though most of the X-ray emis- 
sion was subseque ntly shown to be non- 
thermal in nature. IWinkler fc Long (119971 ) 
obtained a pre-shock density near 1.0 from 
the spatial profile of the X-ray emission, 
but lLong et al.1 tooi ) found n ~ 0.25 
from the value of n e t obtained from the 
Chandra spectrum of their region NW-1, 
which lies immediately behind the region 
of the Balmer line filament we observed. 
The Chandra spectrum showed solar abun- 
dances, which is consistent with the inter- 
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pretation of shocked interstellar gas, and 
the temperature of 0.7 keV is consistent 
with very inefficient thermal equilibration 
between ions and electrons. Thus the 
present results agree well with the results 
of shock wave models of the X-ray spec- 
tra and con firm the pa r amete rs derived. 



However, as lLong et al.l (120031 ) point out 



the shock models did not provide an ac- 
ceptable x-squared fit to the data, so some 
aspect of the physics remains to be under- 
stood. In particular, with n as an inde- 
pendently measured quantity rather than 
a free parameter, reanalysis of the Chan- 
dra spectrum might be able to place better 
limits on the non-thermal emission in the 
NW part of SN1006. 

Another density estimate for the NW re- 



gion o f SN1006 was obtained by lLaming et al. 
(119961 ). who found that the relative inten- 



sities of the UV lines could be explained 
if only about half the O VI emission fell 
within the aperture of the HUT telescope. 
This would require an ionization length for 
the O VI of about 10", or n ~ 0.04 cm" 3 . 
This interpretation is not consistent with 
the values of Uq derived for the same region 
from the HST image, so we conclude that 
th ere is a problem eith er with the models 



of lLaming et al 



(119961 ). or perhaps more 
likely with the reddening correction in the 
far UV. 

A density estimate that is independent 
of SN1006 itsel f comes from H I 2 1 cm 
observations by iDubner et al.l (120021 ). An 
H I feature with a column density of 7 x 
10 20 cm~ 2 at a velocity of -6 km s -1 lies just 
outside the NW rim of SN1006. While the 
velocity does not correspond to the Galac- 
tic rotation value at the distance to SN1006 
(-25 < V LSR < -lGkms" 1 ), theSNR 



is far enough from the plane that the veloc- 
ity could easily differ. IDubner et al.l (120021 ) 
estimate a neutral hydrogen density of 0.5 
cm" 3 , which with the low neu tral fraction 
from Ghavamian et al. rl2002h would im- 
ply a total density of 5 cm -3 , well above 
the values we derive. We conclude that 
either SN1006 is not interacting with the 
cloud identified by Dubner et al., or that 
the shock has only reached the low density 
outskirts of the cloud. 

4.2. Length scales 

The length scale given by the radius of 
curvature is considerably shorter than the 
length of the filament itself. The more or 
less straight portion of the Ha filament 
extends for perhaps 7', or about 4 pc, 
roughly 100 times typical radius of curva- 
ture derived for the ripples. The actual 
line of sight scale of the filament is sev- 
eral times that of the bright rim at the 
trailing edge of the filament. The region 
ahead of the trailing bright rim at position 
28 has a nearly constant surface brightness 
of 2 — 3 x 10~ 5 photons cm -2 s _1 "~ 2 . For 
n = 0.3, f neut = 0.1 and V s = 2900 km s~\ 
this implies an angle of 6 to 10 degrees be- 
tween the shock and the line of sight. The 
roughly constant intensity region is about 
2.8 x 10 17 cm in radial exent, so the range 
of angles implies a depth along the line of 
sight between 1.5 and 2.3 x 10 18 cm, or 
about a tenth the scale of the filament in 
the plane of the sky. The 10 18 cm scale 
corresponds to the smoother curve in the 
schematic diagram in Figure HI while the 
10 17 cm scale derived from the Ha bright- 
ness profiles corresponds to the smaller 
scale ripple superposed to obtain the shock 
geometry shown in the schematic. 
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Even the larger scale inferred for the di- 
rection along the line of sight is several 
times smaller than the length of the fila- 
ment in the plane of the sky. Very faint 
Ha emission, which is not apparent in Fig- 
ure El is seen in the very deep H a image of 
Winkler. Gupta fc Long! (120031 ) extending 
out ahead of the bright filament (their Fig- 
ure 5). This is probably a shock in lower 
density gas that may be farther from tan- 
gency with the line of sight, and it prob- 
ably extends for a distance comparable to 
the length of the filament. 

Figure 2 shows that ripples on the scale 
of 10 18 cm can be seen in the direction 
along the filament, for instance between 
the two boxes where radial profiles were ex- 
tracted. Ripples on the 10 17 cm scale are 
not apparent, perhaps because a 10 17 cm 
ripple with 10% amplitude would be only 
a few resolution elements radially in Figure 
2. However, such small scale rippling un- 
doubtedly contributes to the widths of the 
Ha brightness peaks and prevents them 
from being as sharp as the model peaks, 
as is particularly apparent in the position 
28 profile. It is also quite likely that the 
ripples are not isotropic. If the magnetic 
field lies near the plane of the sky in the 
NE-SW direction, then density structures 
would be elongated in that direction and 
smaller scales would appear along the line 
of sight. There is some evidence that the 
field does lie in this direction based on the 
cap-like morp hology of the non-therm al X- 
ray emission (iWillingale et al.l Il996l ) , but 
there is also a suggestion that the field lies 
in the SE-NW direction based on the rel- 
ative tejrn2eraturesof_protons and oxygen 
ions (IKorreck et al. II2004I ). 

If we attribute the ~ 2 x 10 18 cm 



scale length of the ripples to density fluc- 
tuations in the interstellar gas, we can 
estimate the amplitude of the density 
fluctuation from the amplitude of the 
ripple. The amplitude of the ripple is 
about one tenth the wavelength, and that 
should be about equal to 5V/V. For con- 
stant ram pressure, nV 2 , this requires 
density fluctuations of about 20%. As 
was found for similar ripples in the blast 
wave of the Cygnus Loop, this agrees 
reasonably well with the expectations 
from the spectrum of interstellar turbu 



lence ( iMinter fc Spangled 1 199 71 ; iRaymond 



20031 ). However, the smaller scale ripples 



revealed by the 10 17 cm radius of curva- 
ture also seem to have amplitudes of order 
1/10 the wavelength, and a Kolmogorov 
spectrum of density fluctuations from the 
turbulent cascade woul d lead one to expect 
a sm a ller amplitude (jMinter fc Spangler 



19971 ). iBeresnyak Lazarianl (120061 ) show 



that the spectrum of density fluctuations 
may be considerably flatter than the Kol- 
mogorov spectrum. 

One feature in the filament is probably 
not due to density fluctuations in the ISM, 
however. The bulge near the western edge 
of the filament in Figure [2] is very simi- 
lar to one farther west along the filament 
that coincides closely with a bright knot 
of X-ray emis sion with enhanced elemental 



abun dances (lLong et al.l 120031 ; IVink et al. 



20031). The X- r ay kn ot (position NW- 



2 of Long et all (12003^ is clearly a knot 
of ejecta overtaking the blastwave. The 
bulge at the SW corner of the ACS im- 
age is probably a similar structure with 
lower X-ray contrast due to its smaller size 
and perhaps differences in density, ioniza- 
tion time scale and abundances. According 
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to hydro dynamic simulations of Type la 
SNRs, ejecta knots form at the Rayleigh- 
Taylor unstable contact discontinuity, but 
they re ach at most 87% of the b last wave 
radius (jWang fc Chevalierl 120011 ). Thus 
the ejecta knots should be over a par- 
sec, or about 2' beh i nd th e Balmer fila- 
ment. IWarren et al.l (120051 ) have pointed 



out that if a substantial fraction of the en- 
ergy dissipated by a shock goes into accel- 
erating cosmic rays, the ejecta knots can 
come much closer to the outer shock. In 
the case of the bulge in the southwestern 
corner of the HST image and the X-ray 
knot farther to the SW, the lack of syn- 
chrotron emission at radio and X-ray wave- 
lengths indicates that little energy goes 
into cosmic rays. Instead, since the blast 
wave encounted the dense gas in the NW 
secto r of SN1006 onl y about 180 years 
ago (ILong et al.l 120031 ). the ejecta knots 



have undoubtedly overtaken the blast wave 
because the blast wave decelerated when 
it encountered denser gas. This has an 
important implication for analyses of the 
fraction of shock energy t hat go es into cos- 



mic rays. IWarren et al.l (120051 ) show that 



ejecta knots in Tycho's SNR come much 
clo ser to the blast wave t han predicted 
by IWang fc Chevalierl (120011 ) , and they in- 
terpret this in terms of energy that goes 
into accelerating particles to high energies. 
That argument can be applied to most 
of Tycho's SNR, but regions where the 
SNR shell is flattened or where the Balmer 
line emission is especially bright should be 
avoided in the analysis, because those are 
regions where the shock has probably been 
decelerated by higher density gas. 



4.3. Shock Precursor 

Shock wave precursors have been in- 
ferred from anomalously high line widths 
of the narrow component Ha emission in 
a number of SNRs and attributed to heat- 
ing in the precursor predicted by diffu- 
sive shock acceleration models or to broad 
component hydrogen atoms that over- 
take the shock and heat the upstream gas 
(|Smith et al.lll994j : iHester. Raymond fc Blair 
19941 ) . The small width of the narrow 
component in SN1006 means that there 
is no evidence for such a precursor, but 
it is nevertheless worthwhile to place on 
limit on the emission from such a precur- 
sor. One expects a more or less exponen- 
tial falloff of brightness ahead of the shock 
peak on a scale given by k/V s , where k 
is the comic ray diffusion coefficient, or 
by (noa)^ 1 where a is the charge transfer 
cross section. The former would be about 
1" for k ~ 10 25 cm 2 s _1 . For a speed 
near 3000 km s _1 and the densities in Ta- 
ble 1, the latter would be about 0.6". The 
brightness could be a significant fraction 
of the narrow component brightness, or up 
to half the total brightness of the filament. 

The concave outward ripples we have 
analyzed so far are not appropriate for 
searching for a precursor, because the 
emission from the curved part of the 
shock front could easily resemble emission 
from a precursor. Instead the bright rims 
about 10" ahead are convex outward rip- 
ples, so that any emission ahead of the 
peak would be due to a precursor. Spa- 
tial profiles do show shoulders of order 
l"wide ahead of the leading bright rims 
in several places. However, there is ad- 
ditional faint emission out ahead of the 
main filament along the entire region im- 
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aged flWinkler. Gupta fc Lond (120031 ) . Fig- 
ure 5). That suggests that the leading rim 
might be 5-shaped along the line of sight 
rather than a convex outward arc, so the 
Ha could be similar to that in the models 
described above (Figures 4 and 5). In par- 
ticular, the sections of the filament that 
show an Ha shoulder ahead of the peak 
emission seem to the the regions where the 
diffuse emission ahead of the filament is 
especially bright. Therefore, we can only 
place an upper limit on the precursor just 
ahead of the shock of about 1/3 the peak 
brightness. Stricter limits of about 1/10 
the peak brightness can be placed in some 
sections. Given the lack of evidence for 
heating in a precursor in SN1006, this 
is not surprising. It does suggest, how- 
ever, that morphology alone may not be 
enough to establish the existence of a pre- 
cursor without some additional informa- 
tion such as [N II] or [S I I ] line emission 
(IHester. Raymond &: Blairl |l994j), spatial 
separation of na rrow and broad emission 
( Lee et al. 2006h. supertherma l [N II] line 



widths ( Sollerman et al. 20031 ) 



or narrow 

component line widths of order 40 km s _1 . 

4.4. Bulk velocity contribution to 
line widths 

The widths of the Ha and UV line pro- 
files have been analyzed under the assump- 
tion that the broadening due to bulk mo- 
tions in parts of the filament that are not 
quite tangent to the line of sight is small. 
The range of angles between the line of 
sight and the shock surface derived above 
for the region between the leading and 
trailing rims implies doppler shifts up to 
1/6 the post-shock speed, or V s /8. Mak- 
ing the extreme assumption that both red- 



and blue- shifts are present, and ignor- 
ing the contribution of the bright rims 
that are tangent to the LOS, the contri- 
bution of bulk motions as large as 370 
km s -1 could be present, or a full width 
of 750 km s _1 . With me asured line widths 
of 2290 km s _1 for Ha ( Ghavamian et al. 



20021 ) and 2100 to 2600 km s" 1 for the 



UV lines ( 


Rav 


monc 


. Blair & Lone 


1995; 


Korreck et al. 


2004 


), bulk motions con- 



tribute at most 6% to the measured line 
widths when added in quadrature to the 
thermal widths. 

5. Summary 

An HST image of the Balmer line fila- 
ment in SN1006 resolves the region where 
neutral hydrogen is excited and ionized. 
The thickness of the ionization zone im- 
plies a total density of 0.25 to 0.40cm -3 , 
roughly in the middle of the order of mag- 
nitude range of earlier density estimates. 
It agrees w e ll wit h the density inferred by 
Long et al.l (120031 ) from the ionization time 
scale of the X-ray emission just behind the 
Balmer line filament, providing confirma- 
tion for the shock wave models used to 
interpret the X-ray spectra. Ripples in 
the filament on parsec scales are consistent 
with those expected for a shock propagat- 
ing through a medium with ~ 20% den- 
sity fluctuations on a parsec scale, in ac- 
cordance with measurements of turbulence 
in the ISM. The ripples on 10 17 cm scales 
have higher amplitude than expected for 
a Kolmogorov spectrum of density fluctu- 
ations, however. A bulge in the Ha fil- 
ament corresponding to a knot of ejecta 
shows that ejecta knots can overtake the 
outer SNR shock when the blast wave de- 
celerates upon encountering a dense cloud. 
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The density estimates in this paper 
are based upon model calculations for 
the emission from non-radiative shocks. 
The models took recen t improvements by 
Heng fc McCrayl (120061 ) into account in an 
approximate way, but a more complete 
model would be valuable. Models of the 3D 
morphology of the filament that go beyond 
the simple shape assumed could improve 
the accuracy of the derived parameters, 
especially when spatially resolved spectra 
become available to pin down the addi- 
tional parameters involved. 

This work is based on observations 
made with the NASA HST, and we are 
grateful to Max Mutchler of STScI for 
help in using multidrizzle on the images. 
We thank Dick McCray and the referee, 
Kevin Heng, for very helpful comments, 
and Kevin Heng for computing the aver- 
age velocity of the fast neutrals relative to 
the shock. This work is supported by HST 
Grant GO-10577.01-A to the Smithsonian 
Observatory. This work made use of the 
NASA Astrophysics Data System (ADS). 
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Table 1 

Densities and radii of curvature from Ha spatial profiles 



Position 


RA2000 


^2000 


PA 


n 


FT 


8 


15 ft 02 m 23.488 s 


-41°44'12.45" 


152.15 


>0.35 


<0.4 


9 


15 A 02 m 23.136 s 


_41°44 , 14.62" 


154.65 


>0.40 


<0.4 


10 


15 A 02 m 22.796 s 


-41°44'16.96" 


154.65 


0.30 - 0.35 


0.3 - 0.6 


11 


15 A 02 m 22.450 s 


-41°44'19.29" 


155.90 


>0.35 


<0.5 


26 


15 A 02 m 17.340 s 


_41°44'55.13 // 


153.40 


0.25 - 0.35 


2 - 5 


27 


15^02 m 17.014 s 


-41 44'57.91" 


154.65 


0.30 - 0.40 


0.8 - 1.2 


28 


15 A 02 m 16.682 s 


-41 o 45'00.34" 


153.40 


0.25 - 0.35 


0.7- 1.5 


29 


15 h 02 m 16.340 a 


-41°45'02.74" 


157.15 


0.35 - 0.40 


0.5 - 1 



a 1Q 17 cm 
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Fig. 1. — A three-color figure showing 
HST/ACS F658N (in red) and two energy 
bands of Chandra data (0.4 - 0.7 keV in 
green and 0.8 - 1.5 keV in blue). The field 
of view shown is 279" x 331", with north 
up and east to the left. The X-ray data 
have been binned by 4 pixels and smoothed 
with a 3-pixel Gaussian filter to remove 
pixelation. Note the extended region of X- 
rays coincident with the bulge in the Ha 
image. The drop-off in X-ray at upper left 
is due to a CCD-chip boundary and is not 
real. 

Fig. 2.— The full field of view ACS F658N 
(Ha) image of the NW Balmer filament in 
SN 1006. Two 10" x 20" boxes indicate 
the filament tangencies discussed in detail 
in this paper. 

Fig. 3. — Ha emissivity as a function of po- 
sition behind the shock for a non-radiative 
shock with n = 0.25 cm -3 , f neut = 0.1 
and T e /T p = 0.05 at the shock. This model 
was computed under the assumption that 
charge transfer is isotropi c. To account for 
the ef fects described by iHeng &: McCray 
( 120061 ) we stretched out the broad compo- 
nent emission by a factor of 1.5 when com- 
paring with the observations. The scale as- 
sumes that 1"= 3.1 x 10 16 cm for a distance 
of 2.1 kpc. 

Fig. 4. — Schematic diagram of the ge- 
ometry of the emitting filament along the 
line of sight (LOS). The light line shows a 
large scale ripple in the shock front, and 
the darker curve is the result of adding a 
smaller scale ripple to the light line. The 
fading of the dark curve towards the inside 
of SN1006 indicates the falloff of emissivity 
with distance behind the shock (Figure [3]). 



The trace in the lower right shows the Ha 
brightness as a function of position for lines 
of sight passing near the tangency to the 
LOS. It is obtained by integrating the Ha 
emissivity along each line of sight. 

Fig. 5. — Comparison of observed and 
computed Ha spatial profiles for posi- 
tion 10. The shock wave is concave out- 
ward, with the assumed radius of curvature 
shown in each panel. The inside of the rem- 
nant is to the right. The dashed curves cor- 
respond to the pre-shock densities listed, 
with the lowest density giving the highest 
curve on the inside and the lowest on the 
outside. All the models have been scaled 
to the observed peak intensity by using the 
neutral fraction shown. Neutral fractions 
above 1 are obviously unphysical, implying 
that that combination of pre-shock density 
and radius of curvature cannot account for 
the observed brightness. 

Fig. 6. — Same as Figure 4 for position 28. 
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